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Abstract 

Energy dissipation rate, e, is an important parameter for nearly every experiment on turbulent 
flow. Mathematically precise relationships between e and other measurable statistics for the 
case of anisotropic turbulence are useful to experimentalists. Such relationships are obtained for 
which the measurable statistics are the 3rd-order and 2nd-order velocity structure functions as well 
as the acceleration-velocity structure function. The relationships are derived using the Navier- 
Stokes equation without approximation. Approximate versions are obtained on the basis of local 
stationarity and local homogeneity. The latter are valid for arbitrary Reynolds numbers for the 
case of stationary, homogeneous turbulence. Precise use of the mathematics requires care noted 
in the Discussion section. 
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I. INTRODUCTION 



Recently, many experiments have produced turbulence by novel means and/or mea- 
sured energy dissipation rate e by novel means. Luthi et al. (2005) , use a magnetically driven 
flow that is nearly isotropic, and obtain e from velocity derivatives including 21/ (sijSij), where 
Sij is the rate of strain, 1/ is kinematic viscosity, and repeated indices denote summation over 
coordinate directions. Mann et al. (1999) and Ott & Mann (2005) use an oscillating grid, 
particle tracking, and obtain e from ((a, — a^) {ui — hQ), where Oj and Ui are acceleration 
and velocity. Similarly, Ott & Mann (2000) compare those methods with e determined 
from the 3rd- and 2nd-order structure function. Voth et al. (2002) and La Porta et al. 
(2001) generate turbulence between counter-rotating blades in a cylindrical enclosure; using 
particle tracking, they obtain e from 2nd-order structure functions. Ouellette et al. (2006) 
create the same turbulence and include evaluation of e from ((a^ — a[) {ui — u'^)). Berg et 
al. (2006) use multiple propellers to create turbulence, use particle tracking, and obtain 
e from 2nd-order velocity structure functions. Tsinober et al. (1992), Kolmyansky et al. 
(2001), and Guhtski et al. (2007) measure in the atmospheric surface layer using multi- 
wire probes. Without use of Taylor's frozen flow hypothesis they obtain e from velocity 
derivatives including 2z/ {sijSij). 

A common means of obtaining e is from the inertial range of the energy spectrum or 
that of the 2nd-order velocity structure function. That relationship is based on dimensional 
analysis under the assumption of local isotropy, and based on empirical validation, and 
on empirical evaluation of the Kolmogorov constant. For anisotropic turbulence, that 
empirical basis must be reevaluated for each anisotropic flow. Those relationships have not 
been derived from the Navier-Stokes equation, and are therefore not considered here. 

The usefulness of asymptotic relationships for an experimenter's purpose of quanti- 
fying e using measurements of other statistics is limited. For every different type of flow and 
every Reynolds number, the experimenter must demonstrate to what accuracy the asymp- 
totic hmit permits determination of e. Duchon & Robert (2000) and Eyink (2003) use 
time and space averages of arbitrary extent and v — Q such that the Reynolds number is 
infinite and use an orientation average to remove the effect of anisotropy. Both the time 
and space averages must be of nonvanishing extent (Eyink (2003) and §6.2 in Hill (2006)). 
The requirement v = Q is not applicable to experiments. Eyink (2003) considers the possi- 



2 



bility of experimental tests of their results, concluding that "a slow approach to asymptopia 
makes a direct test of the local results, especially a verification of the numerical prefactor, 
rather more difficult." A likely cause of inaccuracy in experiments of the results of Duchon 
& Robert (2000) and Eyink (2003) is an effect of random sweeping described in §6.1 in Hill 
(2006). 

The time and space averages used by Nie and Tanveer (1999) are not applicable to 
experiments. Nie & Tanveer (1999) base their relationship between e and the 3rd-order 
velocity structure function on time averaging over an infinite duration to remove effects of 
nonstationarity, on space averaging over either all of space without boundary conditions or 
space averaging over an entire spatial period of a spatially periodic fiow to remove effects 
of inhomogeneity, and on a sufficiently large Reynolds number. They use an orientation 
averaging to remove the effect of anisotropy. 

Danaila et al. (2002, 2004, and references therein) present structure function equa- 
tions in which terms that describe particular effects of large-scale inhomogeneity are retained. 
They use data to demonstrate that those retained large-scale terms are the dominant inho- 
mogeneous terms for several types of flows. Their equations are approximate because the 
large-scale inhomogeneity terms are approximated, and the other terms are approximated 
using local isotropy. 

In contrast, the method in this paper avoids approximation. On the basis of algebra, 
calculus, incompressibility, the Navier-Stokes equation, and use of no approximations what- 
soever, the structure function equation that contains energy dissipation rate is given in §4. 
The new, quantiflable definitions of local homogeneity and local stationarity are illustrated 
in §3. In §4, terms that are negligible for local homogeneity and local stationarity are iden- 
tified. For simplicity, but not for necessity, those terms are not carried forward in f|TTl) - fl23l) . 
Those terms are considered in §7. For simplicity, an ensemble average is used herein because 
the slight complications from spatial and temporal averages were thoroughly documented in 
Hill (2002a,b). With attention to those complications, spatial and temporal averages may 
be substituted for the ensemble average used here. An average conditioned on the value 
of some hydrodynamic quantity cannot be used because of its unknown commutation with 
respect to spatial and temporal derivatives; that topic must await future study. 

Several methods of obtaining energy dissipation rate e without use of local isotropy 
are given here. These methods derive from the Navier-Stokes equation with no restriction 
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on the flow symmetry; they differ only in the type of average employed. These methods 
are expected to be of use to experimenters in interpretation of their data. Our purpose 
here is to provide pragmatic means for measuring e using no assumptions and with the 
greatest mathematical precision possible. Cautions to experimenters on precise use of the 
relationships are in §7. 



II. NOTATION 

The velocities, accelerations, and energy dissipation rates at spatial points x and x' 
and times t and t' are denoted by 

Ui = Ui{x, t),u[ = Ui{x', t'), tti = ai(x, t),a[ = ai(x', t'),e = e{x, t),e' = £(x', t'), (1) 

etc.; X, t, x', t' are independent variables. For particle tracking measurements, Ui could be 
the velocity of one particle at position x at time t, and u'^^ could be the velocity of another 
particle at position x' at time t', where t — t' or t ^ t' , but if Ui and u'^^ are the velocities 
of the same particle at different x and x', then clearly t ^ t' . A sequence of point pairs 
may come from particle trajectories, but it is useful to consider them as pairs of coordinate 
locations. Define a new set of independent variables: 



X = (x + x') /2 and r = x - x, and r = |r| ; T = (t + f) /2 and \ = t- t' . (2) 

The significance of variables X and T is that they are the location and time of the measure- 
ment, respectively. Define the following: 

(3) 
(4) 
(5) 

(6) 

Below, we use numerical subscript 1 to denote a component in the direction of r; and 2 and 
3 to denote components transverse to r. 





- u'i) {Uj 






- O 


- U'j) {Uk - Mfe) 


Aij = {di 


- a[) {uj 


- u'j) + [ttj - a'j) {ui - u'. 


F iik = i^i 


- u'i) {Ui 
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III. WHAT ARE LOCAL HOMOGENEITY AND LOCAL STATIONARITY? 



Local homogeneity and local stationarity only apply for very large Reynolds numbers 
and sufficiently small r. Require that the definition of local homogeneity produces the same 
results as homogeneity, and that the approximation be quantifiable. The simplest case is 
the incompressibility relationship on the 2nd-order structure function, namely the divergence 
vanishes: 8^^ (dij) = 0, where we denote an ensemble average by angle brackets. Familiar 
examples are obtained by substitution of the isotropic formula for {dij) into drj (dij) = 0; 
that gives the well-known incompressibility relationship that rdj. (dn) + 2 [(rfn) — (^22)] = 
0; subsequent substitution of the inertial-range 2/3 power-law formula for (dij) gives the 
familiar inertial-range result that (^22) = (4/3) (du), and similarly substituting the viscous- 
range formula for {dij) gives the familiar relationship {{du2/ dxif'') = 2 (^{dui/ dxif''). Use of 
algebra, calculus, and incompressibility, i.e., d^Mi = and d^'Mi = 0, but no approximations 
and no average, gives (Hill 2002a,b) 

dr, [{ui - u'^) {uj - u'j)] = dx, [{ui + u'^) {uj - u'j)] /2. (7) 

Apply an ensemble average and the definition ([3]) in ([7]). Then, to obtain dr^ (dij) = from 
([7]), local homogeneity must be the approximation that 

dr, {d.j) » dx^ {{u, + [u, - u'j) > /2. (8) 

This is a quantifiable approximation because ((uj + u'^) {uj — u'j)) in ^ can be measured 
at several locations X such that the derivative on the right-hand side of ([H]) can be calcu- 
lated numerically. Similarly, local stationarity is the approximation that the derivatives of 
statistics with respect to T are neghgible. The cases of time and space averages applied to 
(IT]) is given in Hill (2002a,b). 

Kolmogorov (1941) introduced a formalism of local homogeneity that uses the 
joint probability distribution function (JPDF) of velocity differences. The moment 
((uj + u'^ {uj — u'j)) in dH]) cannot be calculated from that JPDF. Because of ([7])-(IH]), Kol- 
mogorov's formalism cannot be used to obtain the incompressibility relation drj (dij) = 0, 
nor is it applicable to simplifying the structure-function equations deduced from the Navier- 
Stokes equation (Hill, 2001, 2002a,b, 2006). Kolmogorov's (1941) formalism invokes a region 
of vague size for use of the JPDF. In contrast, ([8]) is truly local because it is a derivative; 
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experimentalists need only displace X sufficiently to use the 3-point numerical derivative 
formula. 

IV. APPLICATION OF THE NAVIER-STOKES EQUATION FOR 
ANISOTROPIC TURBULENCE 

From the Navier-Stokes equations, we obtain an exact equation relating 3rd- and 
2nd-order velocity structure functions and other statistics (Hill, 2002a,b, 2006). "Exact" 
means that no approximations were used; calculus and algebra were used. The trace is 
performed because it greatly simplifies the term involving the pressure-gradient difference. 
We obtain 

An = d-rdii + dx^F uk + dr^duk = 2vdr^dr^dii -2 {£ + £') + W, (9) 

where 

W ^ -2dx^ [{p - p') (u, - u'^] + "^dx.dxAi - ^^dx.dx, (p + p') . (10) 

No average exists in (l9l)- (|T0l) . 

After performing the ensemble average, use of the approximation of local homogene- 
ity as in §3, ie., derivatives with respect to X are negibible, causes the average of all terms 
in (ITOl) as well as the term dx^. {Fuk) to be negligible because they are all the rate of change 
of a statistic with respect to where the measurement is performed, i.e., X. Likewise, the 
approximation of local stationarity in §3 causes the one term (da) to be negligible because 
it is the rate of change of (da) with respect to when the measurement is performed, i.e., T. 
We will return to the evaluation of those neglected terms in the discussion §7. The result 
is the approximate structure-function equation (Hill, 2006). 

(Au) = dr, {duk) = 2udr,dr, {d^) - 2 {e + s') . (11) 

Note the two equality signs in both (Q and (ITT]) . In (fTT]) dr,. {duk) is the divergence of the 
vector {diik); e.g., in Cartesian coordinates {duk) = An {dm) + (9^2 {dii2) + 9^3 {dia). The 
energy dissipation rate in ([9]) and (fTTj) is defined by 

e = 2uSijSij. (12) 

The calculation of exact averages is described elsewhere. (Hill 2001, 2002a,b, 2006) Structure 
functions that contain the two-point pressure difference present formidable experimental 
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difficulties; it is therefore significant that the pressure does not appear in (fTT]) on the basis 
of local homogeneity (Hill 2002a,b, 2006). Local isotropy was not used to obtain (fTTI) ; 
therefore, (II ip can provide methods for measuring the energy dissipation rate in anisotropic 
turbulence. 

For simplicity of notation, let 

e (r) = (5 + e') /2. (13) 

Dependence on r in f|T3|) is because e (r) depends on energy dissipation rates at two points 
separated by r. Anisotropic turbulence causes dependence of e (r) on the direction of r; 
hence the emphasis on r in ( |T3l) . In ( |T3l) the argument list could be (r,^), which assumes 
that the ensemble contains events for fixed r and t. However, if an experimenter chooses 
an ensemble of events having fixed r but various \, then there is an implicit average over 
\ such that t should be deleted from the argument list. Henceforth, t is deleted from the 
argument list. 

For simplicity, consider the case for which r is much greater than dissipation scales; 
for that case we neglect 2h'dr,^dr,, {dij) in (iTTl) . but we later include 2i'dr^dr^ {dij)- Then, 
because A-,- = 2 fa,- — a'A {m — -u'), (fTTI] 



((a, - a'i) {u, - u[)) = -2e (r) . (14) 

For anisotropic turbulence, f|T4l) makes it clear that ((oj — a'^) {ui — u'j)) depends on the 
direction of r. Mann et al. (1999) obtained a relationship similar to f|T^ . and they (Mann 
et al, 1999, Ott & Mann, 2000, 2005) used their relationship to obtain energy dissipation 
rate from their measurements of acceleration and velocity. For our purpose of mathematical 
precision, it is necessary to note the distinctions that (|T4l) is obtained here and in Hill (2006) 
without the assumption by Mann et al. (1999) that a certain derivative moment may be 
neglected and that both e and e' define e (r) in f|T3|) . i.e., two space-time points appear here. 
Average (fT4l) over orientations of r to obtain: 



1 

47r 



((ttj - a[) {Ui - U'i)) dn = -2eorientation (r) (15) 
^orientation (r) = ^ jj sdVt (16) 

which is a function of r; dVt is the differential of solid angle, and the double integral is 
understood to be over An steradians. 
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Use of particle tracking data allows calculation of the average of (fT^ over a sphere 
in r-space to obtain energy dissipation averaged within the sphere as follows: 

\r\<rs 

= -2esphcrc (rs) ■ (17) 

^sphere (ts) =^ jjj Mv. (18) 

|r|<rs 

The sphere has radius rs such that ^sphere depends on rg, not on r. The average over 
orientations of r is performed first, resulting in a function of r; the r-integration is performed 
second. In f|T8|) . the average produces the same result as that given in the first description 
of intermittency theory in Obukhov (1962) which was used by Kolmogorov (1962). 

Another method is to use the equality in ffTTl) that contains dr^. (duk)- To avoid 
substituting data into the divergence of the 3rd-order structure function, the r-sphere average 
is performed and the divergence theorem is used to express the result as an integral over the 
surface of the outward normal of the vector (duk)', then 

3 f frk , , . , 1 




dr, {duk) dr = — ^ / / - {duk) ds=— / / dn 



4:7irgJ J J ' 47rr| J J r ' An 

|r|<rs \r\=rs 

4 

= --rs£sphere (^5) • (19) 

The differential of surface area on the sphere is ds; ^ is the unit vector in the direction 
of r such that ^ {duk) = {dm). Specifically, (dm) = {dm + (^221 + c?33i). The subscript 
notation [{dm)]^^^^^^ means evaluate the quantity within the square brackets at |r| = r^. 
In (1191) . we have the orientation average of (dm). Taylor et al. (2003) use DNS data to 
demonstrate the efficacy of the orientation average of the 3rd-order structure function in its 
relationship to e. 

Neglect of the term vdy^dr, (da) in flTTl) is unnecessary. The Laplacian is the diver- 
gence of the gradient such that the divergence theorem gives 

dvAk {da) jj i^ii) " ^^// [^''1 ^^^^)]|r|=r. d^, 

\r\<rs \r\=rs 
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where dr^ is the gradient in the direction of r. Thus, more generally, for anisotropic 
turbulence 



1 
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{{ai - a[) {ui - u[)) - ydr^dr^ {da) = -2e (r) . 

[{{tti - a-) {Ui - U[)) - ydr^dr^ {da)] dQ = -2£orientation 



^JJ {ia,-a'^iu,-u'^)dn]r^dr 
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(20) 
(21) 

(22) 

(23) 

For stationary, homogeneous turbulence, (|2n]) - (|25]) are valid for all r and all Reynolds num- 
bers. We began with the hydrodynamic definition of energy dissipation rate, e, in f[T^ : 
( IT3l) is only simplified notation. We obtain in ( |20l) - (!23l) three different averages of e. The 
three values differ only because the averaging operations differ. There is no distinction as 
to which value to prefer. Any choice of averaging operation is left to the judgment of the 
experimenter. 



Iff 4 _ 

[{dm) - 2udr^ (4i)]|r|=r-, = "o'^^^sphcrc [rg) 



V. THE LIMIT OF LOCAL ISOTROPY 

We can determine whether or not the above formulas give the corresponding classic 
results for locally isotropic turbulence. All of the energy dissipation rates in ( |20l) - (!23ll 
become the same value for local isotropy, so here we denote them all by (e). For local 
isotropy, both (1201) and (1211) become 

((a, - (m, - u[)) - ur-'dr (r'dr) (du) = -2e, (24) 

where dr^dr,, = r~^dr (r^dr) was used, and fl22l) - fl23l) become 

r 

j {{ai - a[) {u, - u[)) r'dr - ur^dr {du) = -^r^e. (25) 



4 

{dm) - 2iydr {da) = -^re. (26) 

The inertial range formulas obtained from ( |24l) - (!26l) are the well-known formulas 

4 

((oi - a-) (wi - «■)) = -2 (e) and {dm) = --^r {e) 
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The viscous-range formulas from (12^ - (126|) are vr~'^dr (r'^dr) (da) = 2 (e) and udr (da) = 
|r (e). In the viscous range, Taylor series expansion, local isotropy, and incompressibility 
give 

Note that r'^dr (r^dr) = 6. Then, from ([24])-([26D, both ur'^dr (r^dr) {da) = 2 {e) and 
udr (da) = |r {e) give the classic formulas 

VI. DISCUSSION 

The mathematics is precise. Experimenters must be careful to follow with precision 
when using the equations. In particular, precise evaluation of e requires evaluation of 
not only (I20|) - fl23l) . but also the terms in IQ that describe inhomogeneity, i.e., dx,.r uk 
and W, as well as the term that describes nonstationarity, i.e., d'^da. Those terms must 
be operated upon with the same averages that appear in fl2ni) - fl2^ . Danaila et al. (2002, 
2004, and references therein) give approximate evaluations of some of those terms for several 
flows. The pressure that appears in W in ( ITOl) presents a future challenge, but techniques of 
pressure measurement combined with hot-wire velocity measurement have advanced (Tsuji 
et al, 2007). 

The same velocity and acceleration that appear in the Navier-Stokes equation also 
appear above. That is, the velocity fluctuation and acceleration fluctuation do not appear 
above. If the experimenter performs a Reynolds decomposition, e.g. u = (u) + u where 
u is the fluctuation of velocity, then many more terms must appear. Those terms must 
be evaluated quantitatively to obtain quantitative energy dissipation rate. The energy 
dissipation rate may be expressed as being caused by the average flow and the fluctuation flow 
and terms descriptive of the interaction of the two. The subject of Reynolds decomposition 
for the structure-function equations, and the many resultant terms that must be quantified 
or neglected is discussed in Hill (2002b). If those many terms are neglected, the resultant 
approximate equations for structure functions of fluctuations are described fully in Hill 
(2002b). 
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An important advantage of using (An) over dr^ (dak) in (fTT]) is the necessity to 
neglect (9t {da) + dx^ (/^ uk) in © if dr^ {duk) is to be used. For the case of a nonzero mean 
velocity as well as for random sweeping by the large-scale flow, the term dxf, {F uk) can be a 
significant effect [see §7 of Hill (2006) and experimental evaluation in Danaila et al. (2002, 
2004, and references therein)]. 

The definition of local homogeneity used above holds at the center of symmetry of 
a flow. For example, in the center of the cylinder containing the flow between rotating 
blades, the rate of change with respect to where the measurement is performed (i.e., X) is 
zero because any statistic evaluated both above and below the center of symmetry has the 
same value at both points. The same is true for any statistic evaluated at two points equally 
spaced in any direction from the center of symmetry. 

For use of Taylor's hypothesis of frozen flow, the correction for fluctuating convection 
velocity is given in Hill (1996) for any statistic on the basis of local isotropy; that makes that 
correction inapplicable to anisotropic turbulence. However, the correction is qualitatively 
useful. An advantage of the 3rd-order structure function (dijk) is that the correction to 
Taylor's hypothesis caused by fluctuating velocity vanishes in the inertial range of (dijk) 
(Hill, 1996); the same can be shown to be true for (Aij). 

In (!20|) - (pT|) the calculation of the three-dimensional Laplacian of data might pose 
problems. Danaila et al. (2002, 2004, and references therein) have approximated the 
Laplacian by only its r-derivatives on the basis that local isotropy should be adequate at the 
dissipation and viscous-range scales where udr^.dr,, (da) is an important term within ( 120|) - 
( I2TI) . For axisymmetric turbulence one may express the Laplacian in terms of only two 
independent variables using the cylindrical coordinate system used by Lindborg (1995) or 
the coordinates used by Batchelor (1946). The axisymmetric analogues of fl28l) are given 
by George and Hussein (1991). 

The results (l20l)-(!23l) relate the energy dissipation rates on the right-hand sides to 
quantities that are measurable on the left-hand sides. The results are valid for anisotropic 
turbulence. The left-hand sides of ( I20|) -(l23l) are measurable with particle tracking technol- 
ogy (Voth et al. 2002, La Porta et al. 2001, Mann et al. 1999, Ott & Mann 2000, 2005, 
Luthi et al, 2005, Berg et al. 2006, Ouellette et al. 2006). With the exception of the 
acceleration- velocity structure function, the left-hand sides of fl20l) - fl23|) are also measurable 
by multi-wire probes (Tsinober et al. 1992, Kolmyansky et al. 2001, Gulitski et al. 2007). 
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We expect that (l20!) - (l23l) will be of significant use to experimenters. 

The author thanks A. Pumir, B. Luthi, H. Xu & E. Bodenschatz for valuable comments. 
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